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In canonical scalar field inflation, the Starobinsky model (witli a linear potential but discontinuous 

slope) is remarkable in that though slow-roll is violated, both the power-spectrum and bi-spectrum 

. . . can be calculated exactly analytically. The two-point function is characterised by different power 

C^ \ on large and small scales, and a burst of small amplitude superimposed oscillations in between. 

Motivated by string-theory models with stuck branes, we extend this analysis to Dirac Born Infeld 

^^ ' (DBI) inflation, for which generalised slow-roll is violated at the discontinuity and a rapid variation 

^N I in the speed of sound Cg occurs. In an attempt to characterise the effect of non-linear kinetic terms 

5— ( , on the oscillatory features of the primordial power-spectrum, we show that the resulting power 

C*H' spectrum has a shape and features which differ significantly from those of the standard Starobinsky 

■^F ' model. In particular, the power-spectrum now takes very similar scale invariant values on large and 

small scales, while on intermediate scales it is characterised by much larger amplitude and higher 

frequency superimposed oscillations. We also show that calculating non-Gaussianities in this model 

is a complicated but interesting task since all terms in the cubic action now contribute. Investigating 

whether the superimposed oscillations could fit to the Planck Cosmic Microwave Background (CMB) 

data (for instance by explaining the large scale Planck anomalies) with, at the same time, small 

non-Gaussianities remains an intriguing and open possibility. 

^ ; I. INTRODUCTION 

It is by now a near certainty that the universe underwent a period of accelerated expansion — inflation — early 

K>" ' in its history. Indeed, the spectacular Planck data [1, 2] is entirely compatible with Standard Single Field Inflation 

^N ' (SSFI) in the slow- roll regime, and with a canonical kinetic term [3]. However, the nature of the inflaton field still 

2~^ . remains a mystery, and on the theoretical side much work has been carried out in the last years with the aim of 

-^ ' trying to embed the infiationary scenario in a high-energy theory such as string theory. This paper fits into such a 

• , "top-down" approach. In particular our focus is on string motivated models (see below) in which, because generalised 

^^^ ' slow-roll is violated for a few e-folds, oscillations are generated in the power spectrum. Our aim is to determine 

precisely the observational predictions of this class of models. In the process we will characterise how, for a given 

inflaton potential, the properties of the superimposed oscillations in 'P^(fc) — amplitude and frequency for example 

— depend on non-standard kinetic terms in the Lagrangian for the inflaton. 

The model we study is closely related to the so-called Starobinsky model [4. 5] for which the potential V{4)) is linear 
with a sharp change of slope at a certain (J)q : 
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C^ ■ ../ ,s _ j Vo + A+ {(j)- 0o) for (f) > 00, /-, N 

In the following we will take A+ > A_. The change in slope causes a short, of order one in e-folds, period of fast roll, 
and remarkably (in SSFI with standard kinetic terms) both the power spectrum and bispectrum can be determined 
exactly analytically (in all range of parameter space). To our knowledge, this is the only model in SSFI for which any 
exact statements can be made. One finds [4] that there is a sharp rise in the power-spectrum T'f (fc) on scales fc ~ fcg 
(where kg is the mode that left the Hubble radius at = ^o), with 

ita r-Hk)=(^X(^-^)\ lim r-Hk)^(^y(^)\ (2) 
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where Hq is the Hubble scale at the time when (j) = (j)o, and CS denotes the 'canonical' Starobinsky model. Thus the 
increase in power is proportional to AZ — A'^ , and it is followed by small oscillations for fc>fco whose amplitude 
arc rapidly damped out. The different contributions to the bi-spectrum can also be calculated analytically [6, 7] in 
terms of A± and the third parameter of the model Hq, and some ranges of parameter space are ruled out by recent 
constraints on f^-^ from Cosmic Microwave Background (CMB) data [2, 8]. 

Here, motivated by string-theory, we focus on Dirac Born Infeld (DBI) branc inflation (sec e.g. [9-f4]) with action 



5" = 
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X\ 



Ml,R + T(0) Jl + ^ g^'^d^c^dA + V{4>) - T{4>) 



(3) 



where R is the Ricci scalar and M\>\ the reduced Planck mass. As usual r((/)), which we call the 'brane tension' from 
the physical origin of action (3), is determined by the warp factor of the 10 dimensional metric in which the brane 
moves. Provided it is continuous, its precise form is essentially unimportant for this analysis, though to be concrete 
we take 

r(0) = ^ (4) 

which is typical for anti de Sitter warp factors. The potential ^(0) has contributions coming from the interaction 
between the branc and the background, as well as any other branes which may be present in the geometry. Its 
exact expression is not known and here, without attempting to concretely realise the embedding of this model within 
string-theory, we suppose it is given by (1). Such a sudden jump at 00 can be thought of as mimicking the presence 
of a trapped brane stuck at a fixed point of an orbifold symmetry [15], and potentials of this kind have been studied 
in the DBI-literature before e.g. [16]. In particular, if one considers the excitation of particles living on the trapped 
brane, then these can backreact on the inflaton dynamics [17], though the effect has been shown to be small. Note 
that in the limit in which the sound speed Cg — >■ 1, our model reduces to the standard Starobinsky- model discussed 
above. Our model may not be completely realistic from a string theory point of view but its crucial advantage is that 
it allows us to derive explicit analytical results, which can be compared to known results for the standard Starobinsky 
model. We believe that the scenario studied here represents the best compromise between cases in which the string 
model building problem can be properly adressed (often at the expense of solving the equations numerically), and 
over simplified situations in which analytical results can be easily derived. 

Just as in the standard Starobinsky model, in the DBI case, generalised slow-roll is broken for order one efolds 
around 0o • Thus Cg also changes rapidly, and we expect oscillations to be generated in the power-spectrum. How does 
the shape of Vq (k) depend on the non-standard kinetic terms? Are the amplitude and wavelength of the oscillations 
sensitive to the non- linear structure of the Lagrangian? In this paper not only do we determine V(^{k) numerically for 
all Cg , but we also show that the model is essentially completely soluble analytically in the c^ ^ 1 limit, in terms of 
the parameters of the potential A±. Due to the non-linear kinetic terms in the action, V(^{k) differs significantly from 
that of the canonical SSFI Starobinsky model discussed above. Indeed, the action given in (3) now contains a second 
dimensionful potential r((/)) and it is this, rather than A±, which determines the power-spectrum on small and large 
scales: 
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where Tq = T(0o)- As opposed to the canonical Starobinsky model in which there is a sharp rise in power across ko if 
A^ ^ A^ , in the DBI-Starobinsky model there is no rise in power in this case. Rapid variations of Cg occur when the 
field crosses (j)o, and these give rise to large amplitude, high frequency, superimposed oscillations in T'Q{k) which we will 
discuss in section HI. An interesting question is whether or not these oscillations could fit the Planck data, for instance 
explaining the large scale Planck anomalies, while also remaining compatible with constraints on non-Gaussianities. 
Indeed, as we discuss in the conclusions, an interesting new feature of this model is that we expect non-Gaussianities 
to be sourced predominantly from two (or more) coupled vertices, leading to a complicated structure. While it might 
be expected that the constraints are strong at least for Cg ^ 1, for larger Cg the situation regarding non-Gaussianities 
is much less clear while the large amplitude superimposed oscillations remain in the power-spectrum. 

In the context of SSFI, the development of models leading to oscillations was motivated by observed features in the 
power-spectrum of CMB temperature fluctuations [18-22]. Indeed, they may find their origin in initial conditions, 
arising for instance from non-Bunch-Davies initial conditions (see e.g. [19, 23]), or from deviations from slow-roll in 



SSFI. Amongst the models studied in the hterature are, for instance, potentials whose derivative is discontinuous 
[4, 6], as well as potentials which contain a step {e.g. [5, 24]), or a sinusoidal modulation (e.g. [25-27]). Though 
these features can fit data better than a nearly scale invariant power-spectrum (see for example [1, 28]), this is at 
the expense of including extra parameters into the potential meaning that the statistical significance of the features 
is not so obvious [29]. Other than solving for the power-spectrum numerically [30-32], some semi-analytical methods 
have been developed [33, 34] though these are generally valid only in certain limiting cases, for example if the step is 
small or if the scalar field and metric perturbations decouple [27] . 

In models with non-standard kinetic terms, the consequences of rapid variations in Cg in fc-infiation have been 
investigated (both for the power-spectrum and bispectrum) using the effective field theory formalism in [35]. The 
"generalised slow-roll approximation" has been extended to A:- inflation [36], and applied to DBI- inflation in [37], 
though there the authors considered a step-like feature in T{(l>). Notice that Ref. [38] also carries out a detailed 
analysis of the signatures of step-like feature in both T((/)) and V{(j>). 

The structure of this paper is the following. In section II we first discuss the background evolution of the system 
and define the generalised slow- roll parameters (subsection II A). At this point, our discussion is for a general potential 
V{4>). In subsection II B we focus on the Starobinsky model itself, with potentials V{(l)) and T{(j)) given in Eqs. (1) 
and (4) respectively. From the form of V^cj)) we are able to derive exact results regarding the behaviour of the system 
at the transition (j) — (J)q, and these enable us to determine the evolution of all the slow-roll parameters analytically, 
even when slow-roll is violated. Our results, valid in the DBI limit Cg <SC 1, are shown to match perfectly with 
a full numerical solution of the background equations. In section III, using these exact results, we show that the 
calculation the power-spectrum V^^k) reduces to solving single differential equation with particular time-dependent 
coefficients that we specify. The resulting power-spectrum is shown to be in perfect agreement with a full numerical 
determination of the same quantity. We also determine analytically the dependence of 'P^(fc) on A±, on large and 
small scales. Finally we summarise our main results in the section IV, where we discuss in detail the different shape 
of the canonical and DBI Starobinsky power-spectra. We also present a few considerations on the calculation of 
non-Gaussianities in this model and mention interesting directions for future work. 

II. BACKGROUND EQUATIONS AND SLOWLY VARYING PARAMETERS 

A. Exact Equations 

For arbitrary potentials V{(j)) and T((j)), and working in a spatially flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) background geometry with metric ds^ = —dt^ + a^{t)dx^, the Friedmann and scalar field equations of 
motion following from Eq. (3) are given, respectively, by 

Here a dot/subscript (p denotes derivative with respect to cosmic time i/field (j) respectively, and H ~ a/a is the 
Hubble parameter. The Lorentz factor 7 is related to the square-root in (3) and is inversely proportional to the 
sound-speed c^ 

7(0,0)^1= (7) 

Later on we will mainly consider Cg ^ 1, which is the opposite limit to the standard infiationary case, obtained when 
Cg — > 1 (or 7 — > 1). For the moment, however, we leave 7 arbitrary and all the expressions in this section are exact. 

In general, due to their complexity, Eqs. (5) and (6) cannot be integrated exactly unless numerical methods are 
used. However, it is also interesting to have analytical approximations and for this reason, we now define the horizon- 
fiow parameters. While in standard infiation they are defined as the successive derivatives of the Hubble parameter 
[39-41] , in DBI-infiation a second hierarchy of parameters must be introduced in order to describe the evolution of 
the sound speed (or equivalently 7) . This hierarchy is defined as the successive derivatives of the Lorentz factor with 
respect to the number of e- folds N = ln(a/ain) (where a\n is the initial value of the scale factor), see e.g. [42]. Thus 



the slow-roll (or slowly- varying) parameters of DBI-inflation are given by 
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The slow-roll approximation will consist in taking \ei\ <C 1 and \Si\ <C 1 and in sections II B we will see that this 
greatly simplifies the equations describing the evolution of the system. 

For the moment, however, we make no approximation. In order to write down the exact expressions for the first 
few slow-roll parameters, it is useful to note from (5) and (6) that the time derivative of the Hubble parameter is 
given by 
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so that i?? = T((/))(7^ — l)/(4Afp[). From here we can extract a very useful expression for 7, namely 
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which will be used extensively below, and in terms of which derivatives of A^ can easily be calculated 
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The exact expressions for the first flow parameters are then 
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One can also proceed the other way round, and express some important background quantities in terms of the slow-roll 
parameters. For instance, on rewriting T in terms of ei using (13), the Friedmann equation (5) can be re-expressed 
as 

. _ y(0) \^ 27 
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As we will see in the next sub-section, this equation turns out to be crucial in order to understand the behaviour of 
the system. On differentiating with respect to (j) and using the definition of the slow-roll parameters yields 
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A further derivative would give V^^, but the result is somewhat tedious, so we only give it below to leading order in 
slow-roll parameters. Finally, one can also derive a useful relation for the derivative of the brane tension, namely 
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As already mentioned, the above equations are all exact and valid for any potential and any brane tension. In the 
next subsection, we focus on the Starobinsky model itself. 



B. Application to the Starobinsky Model 

We now consider the DBI-Starobinsky model, with potentials V{(f>) and T((/)) given in Eqs. (1) and (4). In this 
case, Eqs. (5) and (6) still cannot be integrated analytically. Nevertheless, as we now show, some exact results about 
the behaviour of the system at the transition can be established. 

1. Exact Results for Starobinsky Potential 

Suppose that the field is rolling down the potential (1), starting at a value cjiin with ^in > 0o- Since V^ is 
discontinuous a.t cj) = 4>o it follows from Eq. (6) that both (j> and (j) (and thus 7) are continuous, but (j> is discontinuous. 
Thus 7 is also discontinuous, and its jump can be read off from the following exact equation [consequence of Eqs. (5) 
and (6)] 



T 
which leads to 



3^70 + 1/^ + (7 -1)^0 , (21) 



[7]± = -f [K^]± = AA^, (22) 

where A^ = A- — A+ < 0. This can be rewritten in terms of derivatives with respect to N as 
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We now study the behaviour of the slow-roll parameters at the transition. From Eq. (13), it follows that ei remains 
continuous. However, since £2 is defined in terms of derivatives of ei which itself contains 7 [see Eq. (13)], the second 
horizon flow parameter €2 is discontinuous. (This is also true in the the SSFI-Starobinsky model.) Moreover, following 
from the definition in terms of 7, the parameter Si is also discontinuous. Its jump across the discontinuity is related 
to that of £2 from Eq. (20): 
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If initial conditions at the beginning of inflation are such that all slow-roll parameters are small, it therefore follows 
that the DBI-Starobinsky model is characterised by a continuous £1 which remains small all the time (hence, inflation 
never comes to an end), and by parameters £2 and i5i which are small far from the transition, but jump and can 
be large at the discontinuity. In this sense, the DBI-Starobinsky model is a direct generalisation of the canonical 
Starobinsky model for which £1 <C 1 and £2 jumps at the transition. The new ingredient is, of course, the presence of 
the parameter Si and we have just seen that this parameter has a jump comparable to that of £2, in particular when 
7 > 1, see Eq. (24). 

2. Integration of Equations of Motion: Numerics and Analytic Approximation 

Having understood the broad behaviour of the background quantities (without using any approximation), we now 
aim to understand their evolution in a more detailed fashion, at the quantitative level. As already mentioned, since 
Eqs. (5) and (6) cannot be integrated exactly analytically, we have to rely either on numerical calculations or on 
approximations. In the following, we use both. 

Let us start with the field 4>. In Fig. 1, we show the exact evolution of the field (we have numerically integrated the 
exact equations of motion) for two sets of parameters and initial conditions. In order to understand this behaviour, 
first notice from Eq. (18) that in the slow-roll regime when £1 ^ 1, and for all 7, the Friedmann equation (18) reduces 
to 
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FIG. 1: Left panel: Evolution of the field for the choice eiin ~ fO"^, e2in ^ 2.5 x fO"^ Sun ^ 1.5 x 10"^, 7in ~ 50, 
(f)in/Mpi ~ 252.6788, Hin/Mpi ~ 3.97 x 10"*^ and No = 10, where iVo is the number of e-folds at which the field goes through 
the transition. This implies Vo/M^i ~ 4.71 x 10"^^ A ~ 6.45 x 10^°, <j)o/Mpi ~ 252.6725 and A+/M^i ~ 1.49 x 10"^^ In 
the following this set of parameters is named "run one". The black solid curve corresponds to A- = 0.01A+, the solid green 
one to j4_ = 0.05^4+, the solid blue one to A- = 0.1^+, the solid pink one to A- = 0.2yl+ and, finally, the solid red one 
to A- — 0.5yl+. The presence of the transition at TVq — 10 is easily visible. We see that the greater the change in the 
slopes, the more the trajectory is modified after the transition. The dark green dashed line represents the slow-roll solution 
given by Eq. (29). The slow-roll trajectory is not affected by the transition because, in the approximation used there, this 
trajectory no longer depends on the parameters A±. If the initial velocity is not very large and the slope change important, 
then the actual trajectory can significantly deviates from the slow-roll, see the example of the black solid curve. Otherwise, the 
agreement is excellent. Right panel: evolution of the scalar field for another set of parameters, namely enn — 7.25 x lO"'^, 
ezin ^ -7.98 X 10-^ Sun ^ 8.27 x 10-^ 7i„ ~ 1723.33, 0in/Mpi ~ 0.7075, -ffin/Mpi ^ 1.82 x 10-^ No ~ 17.24. This implies 
that Vo/M^i ~ 9.98 x 10~^^ A ~ 8.95 x 10^^ <^o/Mpi ^ 0.7070, A+/M^i ~ 4.98 x 10~" and A- = 0.01yl+. In the following, 
we denote this set of parameters by "run two" . The solid black line represents the exact trajectory while the dark green dashed 
line represents the slow-roll solution given by Eq. (29). Despite the fact that A- — 0.01^1+ , as for the solid black line in the 
left panel, the agreement between the numerical and slow-roll solutions, is now very good. This is due to the fact that, for run 
two, the Lorentz factor 7 is larger. As a consequence, the field arrives at the transition with a higher velocity and, therefore, 
is less sensitive to the changes in the slopes. 



It is remarkable that, despite its intrinsic complexity in DBI inflation, the Friedmann equation exactly reduces to its 
standard counterpart when ei <C 1. Thus we also have 
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Now, the dynamics of 4>{N) can be obtained from the exact identity (12) which, for a slow-roll trajectory [using 
Eqs. (25) and (26)], reduces to 
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Notice that this expression is in fact valid for any potential and any brane tension provided the slow-roll approximation 
holds. It was established for the first time in Ref. [43]. 

For the DBI-Starobsinky model with potential given in (1), Eq. (27) yields 
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where we have chosen the minus sign since for the Starobinsky potential N increases as the field rolls down the 
potential towards smaller 0. Notice that the origin of the square-root is the 7 factor in (12). The above expression 



is still too complicated to allow an exact integration to determine the field trajectory: further assumptions must be 
made, and the first we make is to assume vacuum domination Vq ^ ^±("0 ~ 4'o) for ^-U ''!'■ Then the trajectory is 
given by 
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which results in Elliptic functions. Since this is still not especially illuminating, we make a second assumption, namely 
that we work in the DBI regime 7 ^ 1, or equivalently Cg ^ 1. Referring to (11) this implies that we neglect the "1" 
in 7, which in Eq. (28) translates into 
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The solution is 
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where A^o denotes the number of e- folds when (j> reaches the transition at (pg. The shortcoming of this expression is 
that, because of our successive approximations — vacuum domination, slow- roll ei <C 1, and DBI-regime 7^1 — we 
have lost the dependence on the coefficients A± . The exact evolution of the scalar field is compared to the slow-roll 
trajectory (29) in Fig. 1. 

We now study Lorentz factor 7 given in Eq. (11) in more detail. Its exact (numerically solved) evolution is shown 
in Fig. 2. If the slow-roll approximation is satisfied (that is to say far from the point where the derivative of the 
potential is discontinuous), then 
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where the last expression is valid in the vacuum dominated regime. Furthermore, in the DBI regime 7 3> 1, the 
second term in the square root must dominate so that 
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where </>(A'') is given in (29). Notice that the Lorentz factor does depend on A±. 

During the transition, the above expression is clearly no longer valid since slow-roll is violated. But we can determine 
the evolution of 7 through the transition using Eq. (23) which, in the 7^1 limit, reduces to 
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Therefore, during the transition era, the Lorentz factor decreases exponentially (recall that /S.A < 0) and 
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since the continuity of 7 requires that 7'^'(A'o) = 71r''(^o)j while 7(A^ ^ A^o) = TsrH^)- Once again, this excellent 
fit for 7 is shown in Fig. 2. 

We now turn to the behaviour of the slow-roll parameters. The quantity ei is determined directly from Eq. (13). 
Away from the transition, where the slow-roll approximation is valid, one obtains 
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while during the transition Eq. (13) yields 
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FIG. 2: Top left panel: Evolution of the Lorentz factor 7(A'') for "run one" (see the definition of "run one" in the caption of 
Fig. 1) for A- = 0.01yl+ (solid black fine), A- = 0.05A+ (sohd green line), A- = 0.1A+ (solid blue line), A^ = 0.2A+ (solid 
pink line) and A- = 0.5^+ (solid red line). The transition at A^o = 10 due to the discontinuity in the slope of the potential 
is clearly visible in this figure. Top right panel: "run one" for A- = 0.01A+ (same as the sohd black line in top left panel). 
The dashed green curve corresponds to the approximate slow-roll evolution of the Lorentz factor given by Eq. (30) and valid 
after the transition. The dotted red cu rve corre sponds to neglecting the term one inside the square in this expression and, 
therefore, to assuming that 7sr(-'V) = y''X/{3Vo)A-Mpi/(j>^{N). As shown in the inset, the dashed green line is an excellent fit 
while the dotted red line is not accurate enough. This is because the case A- — 0.01^4+ corresponds to a brutal change in the 
slope of the potential such that the field velocity strongly decreases after the transition. As a consequence, the Lorentz factor 
approaches one and the factor one in the square root in Eq. (30) can no longer be neglected. Bottom left panel: same as 
top right panel but with A- — 0.2^+. As shown in the inset, this time, both the dashed green line and the dotted red line are 
good fits of the numerical solution. Clearly, this is because the change of slopes is less abrupt and, therefore, the field velocity 
decreases less at the transition. As a consequence, the Lorentz factor remains large compared to one and the factor one in the 
in the square root in Eq. (30) can now be safely neglected. Bottoni right panel: same as top right panel but for "run two" 
(see the definition of "run two" in the caption of Fig. 1). The dashed green line and dotted red line are excellent fit of the 
actual numerical solution (see the inset) despite the fact that the change in the slopes is abrupt, A- = 0.01 A+. The reason for 
this behaviour is of course that the initial value of the Lorentz factor is higher. 
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FIG. 3: Left panel: Evolution of the first slow-roll parameter ei for "run one" (see the definition of "run one" in the caption 
of Fig. 1) for A- = 0.01A+ (solid black line), A_ = 0.05yl+ (sofid green line), A- = 0.1^+ (solid blue line), A- = Q.2A+ (solid 
pink line) and A- — 0.5^4+ (solid red line). Right panel: Evolution of the first slow-roll parameter ei for "run two" (see 
the definition of "run two" in the caption of Fig. 1). The dark green dashed line corresponds to the "slow-roll" solution given 
by Eq. (35) (and valid only after the transition). As can be observed from the figure, this is an excellent fit to the numerical 
solution. 



Thus ei also decays exponentially during the transition, from a very small value to another small value: in other words 
ei is small even in the transition region around (Jjq. These analytical estimates are compared to the exact evolution 
of ei in Fig. 3 and one again notices that the matching is excellent. 

Let us now study the slow-roll parameters 6i and 62. Upon using Eq. (15), one obtains the following expression 
valid only far from the transition 



As opposed to ei , this slow-roll parameter does not depend on A± . During the transition, one has 

3AA e-3(^-^o) 



s[-\n) = 



A_ l-{AA/A-)& 



-3{N-No) 



SlsniN). 



(36) 



(37) 



The behaviour of the slow-roll parameter 62 can be obtained in the following way. Starting from Eq. (26) it follows 
that 






7 



(4ei - 62 - (5i) 



Thus in the slow-roll regimes on cither side of 0o where V^^ = 0, 



4ei(±) 



£2 



(±) 



■-^i. 



(38) 



(39) 



On the other hand, we have established that during the transition, £2 
combining Eqs. (37) and (39), we find 



.(-) 



{N) 



iAA 



-3(Af-Wo) 



A- l-(A.4/A_)e-3(^-^o) 



(5i when 7 > 1 [see Eq. (24)]. Thus on 
+ Ae,\rJ{N)^5,,,{N), (40) 



where ei'~^ and 5i^^ are given in Eqs. (34) and (36) respectively. The above considerations are checked in Fig. 4 
where the above analytical estimates are shown to be excellent approximations to the exact numerical evolutions of 
the two slow- roll parameters £2 and 5i. 



10 





_ 


1 1 1 





. 


: 










y"^ /^/^ / j^ 






/ / y / / 






/ / r ^/ f / 


. 










- 


1: \ 


1 


1 


1 








////' / -> 




\ 


'i 


- 


/ / / / / 


1 


- 






/// / / _g 


f .; 








1 












\ y ^ ^^.---^ 


" 






;__,>^^ 1 5 10 15 20 25: 


•i 


- 








1 1 1 ■ 





-1 


■ 1 


: 


'- < 




-2 
-3 


/ a, ■ 

J 6, analytical ■ 

' 



10 11 12 13 

e-folds number N 



14 



16 



17 



18 19 20 
e-folds number N 



21 



22 



FIG. 4: Left panel: evolution of the slow-roll parameters £2 and (5i for "run one" (see the definition of "run one" in the 
caption of Fig. 1) for A- — 0.01A+ (solid black line for £2 and dashed black line for 5i), A^ — 0.05 A+ (solid green line for £2 
and dashed green line for 5i), A- = 0.1^+ (solid blue line for £2 and dashed blue line for Si), A- — 0.2^1+ (solid pink line for 
£2 and dashed pink line for 5i) and A- = 0.5^1+ (solid red line for £2 and dashed red line for 5i). Outside the region of the 
transition, around A^o ~ 10, the evolution is featureless as shown in the inset. When the change in the slope is not too abrupt, 
one observes that £2 ~ 5i . Of course the most important result shown in the figure is that £2 and 5i can be large during the 
transition. When the kinetic term is standard, slow-roll violation corresponds to a situation where £1 ^ 1 (i.e. infiation never 
stops) and |£2| > 1. In the DBI case, this situation generalises to £1 <^ 1, \e2\ > 1 and \Si\ > 1. Right panel: evolution of 
the slow-roll parameters £2 and Si for "run two" (see the definition of "run two" in the caption of Fig. 1). In this case we 
have £2 — Si with a very good approximation despite the fact that A- = 0.01^4+. As before, this is due to the fact that the 
initial velocity of the field is much larger than in the left panel. The dark green dots correspond to the "slow-roll" solution 
given by Eq. (37) (and valid only after the transition). As can be noticed in the figure, this is an excellent fit to the numerical 
solution. The small peaks around A^o are simply numerical artifacts whose origin stems from the fact that we have modelled 
the Heaviside function with a hyperbolic tangent. The amplitude of these peaks can be decreased at will by increasing the 
sharpness of the hyperbolic tangent. 



Finally, we study the behaviour of the quantities 6162 and 6263. These two combinations are important because they 
appear in the effective potential for the cosmological perturbations (see the next section). Just after the discontinuity, 
where €2 — Si ^ ei, we find from Eqs. (16) and (17) that 



6162 
£263 



-3^1 - Sf 
-3e2 - 4: 



T00 
m 



-362 



(41) 
(42) 



where we have neglected the T^^ term in S1S2 which is small. Thus, substituting (40) yields the fit after the disconti- 
nuity of 



6162 



£263 



-9- 



AA 



.,-3{N-No) 



A^ 



"l-(A^/A_)e-3(^-^o)l 



(43) 



This analytical approximation is tested in Fig. 5. As can be seen in this plot, the matching to the exact numerical 
solution is excellent. 

To summarise, the above considerations show that on assuming that potential is vacuum dominated, that the system 
is initially in slow-roll, and that 7^1 for all times, then it is possible to obtain excellent analytical approximations for 
each relevant background quantities and slow-roll parameters throughout the evolution, despite the fact that slow-roll 
is violated at the transition. This is important since the slow-roll parameters control the evolution of the time- 
dependent frequency of each Fourier mode of the perturbations. In the next section we study in detail cosmological 
fiuctuations, our final goal being to determine their power-spectrum. Before doing so, however, we end by noting 
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FIG. 5: Left panel: evolution of SiS^ for "run one" (see the definition of "run one" in the caption of Fig. 1) for A- = 0.01A+ 
(solid black line), A- = 0.05A+ (solid green line), A- = 0.1 A+ (solid blue line). A- = 0.2 A+ (solid pink line) and A^ = 0.5yl+ 
(solid red line). Right panel: evolution of S1S2 (solid black line) and £263 (dashed red line) for "run two" (see the definition of 
"run two" in the caption of Fig. 1). We observe that S1S2 — £2^3. The dark green dots correspond to the "slow-roll" solution 
given by Eq. (43) (and valid only after the transition). As can be noticed in the figure, this is an excellent fit to the numerical 
solution. The vertical lines at A'o originate from the Dirac function in V^^ in Eq. (16). 

that in the slow-roll regime far from the discontinuity, 7^^^ as well as all the slow-roll parameters, only depend very 
weakly on N: indeed they remain essentially constant on either side of the discontinuity (as can be seen in Figs. 1 - 
5). The origin of this behaviour is the small variation of cj) with TV as seen in Fig. 1. For this reason, in the following, 
we will often omit to write the explicit N dependence of the slow-roll parameters in the slow-roll regime. In the above 
expressions, this amounts to approximating (j){N) by 0(A^o) = (f'Oj so that for example from (34) 



ei 



(±) 



MpiA± / 3 



while from Eqs. (31) and (33) 



7. 



(±) 



-^^ i'-HN)^,i-j 



3Vr 



AA 
Al' 



'3{N-No) 



1 



(44) 



(45) 



Having completed the study of the background, wc now turn to the perturbations, in particular to the calculation of 
the two-point correlation function. 

III. POWER-SPECTRUM 

In this section, we are interested in scalar perturbations. It is well-known that they can be characterised by a single 
variable, the so-called Mukhanov- Sasaki quantity, v{rj,x). Its Fourier amplitude obeys the equation of a parametric 
oscillator, namely 



— ]Vk 

7 



0, 



(46) 



where the prime now denotes a derivative with respect to conformal time, k denotes the comoving wave number of 
the Fourier mode under consideration, and z is given by [43] 



2(77) = a(?7)Mp 



(47) 
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The effective potential z" j z which determines the evolution of the scalar perturbations can be determined directly 
from Eq. (47) and the definition of the slow-roll parameters in Eqs. (8) and (9): its exact expression is 



z 



= a^H^ 



2 - ei + -62 + -62 - 2*=i^2 + 2^2e3 + (3 - ei + 62) 5i + 5l + (5i(52 



(48) 



(Note that if 5i ^ 82 ^ 0, this reduces to the SSFI expression, see eg. [6]). Another important difference with respect 
to conventional inflationary theory is the presence of the term 1/7^ = c^ in front of fc^ in Eq. (46), responsible for 
the fact that perturbations now propagate with a speed c^ different to the speed of light. As a result, the question of 
how initial conditions are chosen is more subtle — indeed, one usually assumes that the initial state is the adiabatic 
vacuum for which 



"'"'"^^^JfM)"- 



±i I a;(/c,T)dr 



(49) 



''{k,r]) = c^k'^ — z" jz. The adiabatic approximation is vahd if IQ/w^j <?; 1 where Q = 3a;''^/(4w'^) —u!"/{2lo). 
In the standard inflationary context, the modes are initially within the Hubble radius and u) c^ k. It is then obvious 
that the Wentzel-Kramers-Brillouin (WKB) approximation is valid. In the DBI case, however, the modes arc initially 
within the sonic scale and w ~ Cgfc. Because of the non-trivial time-dependence of Cg, it is not obvious that jQ/w^l ^ 1. 
If this is the case, then one simply looses the ability to choose a well-defined and well-motivated initial state. In fact, 
it is easy to show that, on sub-sonic scales, 

"'^' ^<5i - ei<5i + 5,52 + 151] . (50) 



We have seen that, in the DBI-Starobinsky model, all the slow-roll parameters are small initially {i.e. far from the 
transition). It follows that, in this model, one can identify a well-defined initial state Vk{ri) = 7/(2fc) e^''^''/'''. This 
choice is made in the remainder of this article. 

The power-spectrum, or two-point correlation function, is defined by 

where C,k = Vk/ z is the curvature perturbation, and the right hand side is evaluated in the limit, kc^/{aH) ~ —krj/^ — >■ 
0. To find P^{k), we must integrate (46) for each mode starting from the initial conditions discussed above. In general 
this cannot be done analytically due to the complexity of the equations. As a result, one possibility is to integrate the 
system numerically, and for this purpose we have written a numerical code which exactly integrates the background 
and the perturbations mode by mode (and for arbitrary values of the parameters of the model, and hence any Cg). The 
result is displayed in Fig. 6 (solid blue line). We see that, as k increases, the power-spectrum rapidly dips and reaches 
its minimum at a scale fc* which is a large multiple of kg = —l/rjQ (namely the mode which left the Hubble radius at 
(f> ~ ^o)-^ This dip is followed by large amplitude, high frequency oscillations. Notice that for the parameters chosen 
in Fig. 6, the power-spectrum takes the same scale-invariant value on large and small scales. In the conclusion, we will 
compare this power-spectrum with that obtained in the canonical Starobinsky model for precisely the same values of 
A±: as we will see, although the shapes share the same general aspect, they are in fact very different. 

It is also interesting to have an analytical expression for the power spectrum since this can help understand how its 
shape is modified when the parameters of the model are changed. Therefore we now aim to develop approximation 
methods in order to address this question. To solve Eq. (46) and hence determine the v^ and Vc^{k), we use the 
results of section II B in which the slow-roll parameters were determined as a function of N , and hence as a function 
of conformal time 77 since 77 = rjQe~^^~^°' . Our analytical approximation will therefore only be valid in the Cg ^ 1 
regime for which the results of section II B hold. Before the transition 77 < 770 , the slow- roll parameters are small and 
Eq. (48) reduces to z" / z ~ 2a'^H'^ ~ 2/77^ so that the solution of (46) is 



^•"^Vf^-f )'""""■"■ <^^' 



^ Notice that wc have chosen to define fco without any factor of 7 despite the fact that, physically, the relevant scale is aH/c^ . The reason 
for this is that, in this model, 7 varies significantly with time and it is not very convenient to include it in the definition of the prefered 
scale. 
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FIG. 6: Power spectrum of the DBI-Starobinsky model and various approximations. The sohd blue line represents the exact 
power spectrum obtained with a numerical, mode by mode, integration for "run two" . The dashed red line gives Vc^ obtained 
from an integration of Eq. (55). The dotted black line corresponds to the power spectrum derived from Eq. (62) and is valid 
only in the large scale limit while the dotted dashed green line comes from Eq. (60) and is relevant only in the small scale limit. 
Finally, the grey horizontal line gives the asymptotic values of the power spectrum, see Eqs. (66) and (64). 



where 7^+-' is given in Eq. (45). After the transition, the effective potential z" / z can be calculated using the fact that 
(5i ~ 62 ^ ei while 5i52 and 6162 are given in Eq. (41) and (42) respectively. On substituting, one finds that the terms 
linear in slow-roll parameters in (48) cancel (just as in the canonical Starobinsky model), but the quadratic terms do 
not, resulting in z" / z ~ 'h?{2 + ie^/A). Thus after the transition, the mode function satisfies 



1 



72(77) if 



\4iri) 



= 0. 



where, from Eq. (40) 



£2 



.(-) 



3AA 



-Z{N-No) 



= 3- 



Lorj'^ 



A^ l-(A.4/A„)e-3(^-^") (1- (1.7/3) 



with uj = {/S.A/A^)/rjQ > 0. Thus, on using Eq. (45), we finally arrive at 



vZ" + 



27 



~ 2 4 



7^hM {I -Coif) 



r,,i3^2 



77^ 



4 (1 - 0)7/3) 



3^2 



0. 



(53) 



(54) 



(55) 



This equation is one of the central results of this paper, since the power-spectrum is determined directly from its 
solution. In order to find ■p^(fc), the modes v^, and their derivatives must be matched at 77 = 770 = — l/(ao-ffo) [so 
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FIG. 7: Approximations for an analytic solution of the mode equation (55). Solid lines: approximations for 7 (red) and £2 
(blue) made on small scales. On large scales, the 7)-dependence of £2 is taken into account (dashed line), namely from Eq. (40), 



£2 



{3AA/A^)e-'"-'^-^°^/ [1 - (Ayl/yl_)e"^(^-^«)l. 



as to determine the two integration constants associated with Eq. (55)]. This can be done by carefully considering 
of the behaviour of z and z" /z across the transition. Indeed, on recalling that z = Mpia{ri) \/1e\ (j])^{r\) , and using 
Eqs. (33) and (35), one finds 



z 2 A+ 770 2 A+ 



Hence the matching conditions at the transition are 



"fc M = «fc ('?o) , < iva) - v'kim) = o-aHo^—^Vf^ (ryo)- 



9AA 

•2Z 



(56) 



(57) 



Unfortunately, Eq. (55) is not soluble analytically. However, it can be solved numerically. At this point, one could 
wonder whether we have gained something given that our aim was to derive approximate analytical formulae and 
that we have already determined 7^^ exactly by means of a mode by mode integration. However, integrating a single 
differential equation is much easier than writing a mode by mode numerical code and, moreover, as we will discuss 
below, Eq. (55) can also be approximated analytically. In Fig. 6 we plot the power-spectrum obtained by solving 
(55) numerically, and then substituting into (51) (red dashed curve). This is compared with the fully mode by 
mode numerical calculation of the power-spectrum (solid blue line) obtained by solving both the background and 
perturbation equations numerically. The agreement is excellent, thus validating all the reasoning used to arrive at 
(55). 

Our aim is now to determine analytically the dependence of the features described above on the parameters of the 
model A±. The correct behaviour of 7^^(fc) on small and large scales can be captured analytically via the following 
approximation scheme. First introduce a conformal time rji > rjQ so that the period after the transition is split into 
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two parts, see Fig. 7. When t] > iji, which we denote by 'region 
Eq. (53) in is then particularly straightforward, namely 



we take e2 ~ and 7 



'SR 




1-i 




7'-" 

,• 'SR 



1 + i2^ ) e^'^^/'^sR* 
kr/ / 



The solution of 



(58) 



The integration constants a\^ and /3k, to be calculated below, determine the power-spectrum since from Eq. (51) it 
follows that 



Vdk) 



i?2 7(-) 

— 'SR 

87r2M2j £- 



|/3fe -Ckfcl 



27r 



Hi 

Tn 



\Pk 



ak\ 



(59) 



where Tq = T{4)q). In the intermediate region, "region I" for which rjQ < rj < iji, we make the following assumptions 
(see Fig. 7). On small scales, we assume £2 constant, with 62 — e2(??o) = 3AA/A+ as obtained from Eq. (40), while 
from Eq. (33) the evolution of 7(77) is approximated by 7 ~ 7^^' iv/vo) ■ Continuity of 7 at 771 determines 771 to be 
given by 771 = ^ITsrVTsr')^^"^- ^^ ^^^^ case, Eq. (53) has the exact solution 

^.small = \[^ \C^{k)Mx) + C2{k)Y,{x)\ (60) 

where Ci{k),C2{k) are /c-dependent integration constants, and 



x{r),k) 





(61) 



In the intermediate region and on large scales, our assumptions are slightly different. Since 7^1, here we expect 
^2/7^ ^ 1, and hence we simply first neglect this term Eq. (55), while keeping the exact 77-dependence of £2. We can 
then integrate exactly to find 



V 



+ Ci{k)r^'\-2 + Cor]"^) 



(62) 



where C3{k) and C4{k) are again fc-dependent integration constants. The different integration constants Ci{k), 
I = 1, . . . , 4 are determined from the boundary conditions at 770 which, from Eq. (57), read 



^('yo) = v^irjo), vl'{r]o)~v+'{7]o) =aoHo-—v+{7jo), 



(63) 



where the mode solutions in the "+" region, v'j^ , are given in Eq. (52). In turn, the integration constants af, and 
/3fc are then obtained from matching at 771, where the relevant conditions are vj,{rii) ~ f^ (771) and ^^'(771) = ^^^'(771), 
namely the mode functions and their derivatives are continuous. Finally, the power-spectrum 'P^(fc) is then obtained 
from Eq. (59). 

Carrying out this process on large scales using (62) gives an analytic expression for V({k), which is plotted in Fig. 6 
(dotted black line). As can be seen, this is good fit to the full numerical solution (solid blue line) on large scales. 
Since the analytical expression is somewhat long, and so we do not give it explicitly — however, we can use it to 
extract the dependence of V(^{k) on the parameters of the model. In particular, we find that 



fc/fco^o ''^ ^ \27tJ \To 



1- 



2A] 



A] 



(64) 



which agrees very well with the numerical values found for different A±. Notice that when A-/A-^- -^ 0, on large 
scales the power-spectrum tends to the asymptotic scale invariant value {Hq/2tt) (^Hq/Tq), drawn by a grey dotted- 
dotted-dashed line in Fig. 6. With our analytic approximation for Vi^{k), one can also determine the position of the 
first feature in the power-spectrum, namely the smallest value of fc = fc» for which V(^{k^) — > 0. This reduces to 
solving an implicit equation of the form tan {k^Tji / j^~>) = f (fc*77i/7^^') where the function f{x) is a ratio of two 
polynomials, determined straightforwardly from the matching, and with A± dependent coefficients. One can check 
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that the analytic result agrees well with that obtained from a numerical solution of the full power-spectrum. It can 
be approximated by 



ko 



A- 



0.81 



(65) 



thus showing the dependence on A±: the larger A^/A^ the more the feature shifts to large scales. 

Once again, carrying out the same procedure on small scales using (60), we again determine an analytic expression 
for Vc;{k) now containing Bessel functions (which are responsible for the oscillations in the power-spectrum). Again 
the expression is rather complicated, but it is a good fit to the numerical result (see Fig 6, dotted dashed green line). 
We can use it to extract the fc — >■ oo behaviour of the power-spectrum (this limit is particularly simple since then the 
arguments of the Bessel functions are large, and they can be replaced by their appropriate asymptotic expansions). 
We find that the scale invariant value of the power-spectrum on small scales is given by 



hm Vcik)=(^X(^ 
fco-i-oo V 27r / \ To 



fc/fc, 

as advertised in the introduction, and also shown in Fig. 6. We can go further: when A_ <^ A^ it follows that v : 
and then the exact expressions for the Bessel functions are known. Our analytic result then gives 



p(^"«^+)(fc) = 




sin[26'(fc)] 
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(for /O fco) 



(66) 
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(67) 



where 
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for 7'+^ > 7'"). 

'SR ^^ 'SR 



(68) 



2/3 



Thus the wavelength of the oscillations is approximately given by ^ 7^+) (^srVTsr^)" " ^^'^ depends on A± through 
the dependence of 7^^-' on these parameters. Relative to the asymptotic value of V({k — )■ 00), the amplitude of these 
oscillations is determined by the ratio 37^"'"'fco/^- When 37'"'"'fco/'>^ > Ij it is quadratic in this parameter, while when 
•^^sR^ka/k < 1, it is linear. This can be observed in Fig. 6 where, for fc/fco — 37^^-' ~ 5 x 10^, the slope of the envelope 
can be seen to change. 



IV. DISCUSSION AND CONCLUSIONS 



The main purpose of this paper was to study the signatures of the inflationary DBI-Starobinsky model, for which 
the potential is linear with a sharp change in slope at a certain (j)Q and the kinetic term a non-minimal, DBI, one. 
In the case of canonical inflation with such a potential, both the power-spectrum V(^{k) as well as the bi-spectrum 
are exactly soluble analytically [6, 7]. Here we have addressed the following questions: what is the shape of V({k) in 
the DBI-Starobinsky case? What signature do the non- linear kinetic terms leave? Does V^ik) still rise sharply from 
small to large scales? 

To approach this problem, first we studied the homogeneous background evolution of the field and generalised 
slow- roll parameters. We showed in section II B that in the DBI regime 7^1, these quantities can all be determined 
analytically to very high accuracy. We also showed that this model is characterised by a first slow-roll parameter 
ei which is tiny throughout the evolution of the system, while the other generalised slow-roll parameters all become 
large after the transition at (/)o, decaying back to small values over a few e- folds. Armed with the analytical expres- 
sions for the slow-roll parameters, we showed that the power-spectrum can simply be obtained by solving the mode 
equation (55). Furthermore, a numerical solution of this equation was shown to agree exactly with a fully numerical 
determination of the power-spectrum of the model (both at the back-ground and perturbative level), see Fig. 6. 

It is also interesting to compare our results to the predictions of the canonical Starobinsky model, see Fig. 8. 
Following an analytical approximation of Eq. (55), we were able to show that in the DBI Starobinsky model, it is 
actually the second dimensionful potential T[<p) (rather than A±) which determines the power-spectrum on small 
and large scales, as summarised in Eqs. (64) and (66). Thus as opposed to the canonical Starobinsky model in which 
there is a sharp rise in power across fco if A_ ^ A^^ in the DBI-Starobinsky model there is no rise in power in 
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FIG. 8: Comparison of the DBI-Starobinsky and CS-Starobinsky power spectra for the same values of A+ and A^ corresponding 
to "run two" . Obviously, the amplitude of the oscillations is much smaller in the CS-Starobinsky case but there is a rise in 
power that does not exist in the DBI-Starobinsky power spectrum. 



this case, see Fig. 8. We were also able to determine the position of the first feature in the power-spectrum, given 
to a good approximation by k^,/ko ~ {A+/A-)p with p ~ 0.8. Again this differs from the CS model in which this 
feature is at fc* ~ fco: independently of A±. Finally, in the A_ ^ A^ limit, we have shown that the wavelength of 
oscillations in the power-spectrum does depend on A±, as opposed to the CS model, see Fig. 8. Furthermore, relative 
to the asymptotic value of V(^{k — > oo), the amplitude of oscillations is now much larger — rather than being of order 
|AA/^+| -- 1 (in the CS model), it is now of order 7^+^ > 1. 

The next step would obviously be to compare in detail the DBI-Starobinsky model to CMB data, and particularly 
the recently released Planck data. This would require interfacing the numerical code used in this paper to calculate 
the power spectrum to a CMB code (typically the CAMB code [44]), and then in turn to a code allowing us to explore 
the corresponding parameter space (typically the CDSMDMC code [45]). Moreover, we would also need to include in the 
analysis the constraints coming from the higher correlation functions (see below). Clearly, this is beyond the scope of 
the present article since, here, we mainly focus on the physical properties of the system rather than on data analysis. 
It is interesting, however, to have a broad idea about the physical values of the parameters. To this aim, we have 
represented in Fig. 9 the Planck "step model" best fit (solid blue line) [1] with the DBI-Starobinsky power spectrum 
(solid red line) with 7in ~ 5 and A_ ~ 0.92A+. Notice that we are now, therefore, in a very different regime to that 
discussed above since A^ is similar to A- and 7 is not very large. In this plot, the overall normalisation is arbitrary 
(we have normalised the two spectra differently on purpose in order to make easier the comparison of the two shapes) . 
Notice that the x-axis is now the physical k today and, hence, the position of the first feature is essentially arbitrary. 
In particular, it will depend on the post-inflationary evolution. On the other hand. Fig. 9 also indicates that the 
shape of the DBI-power spectrum is not the same as the CS-Starobinsky model and, as a consequence, a rigorous 
Bayesian exploration of the parameter space of this model seems to be required before one can conclude whether the 
model is ruled out or, on the contrary, whether it could explain the Planck anomalies and, therefore, improve the fit. 
In any case, from the above discussion, it is clear that the ratio A^/A_ cannot be too different from one otherwise the 
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FIG. 9: Comparison of the DBI-Starobinsky power spectrum (solid red line) with the Planck "step-inflation" model best fit (solid 
blue line). It should be noticed that the DBI-Starobinsky power spectrum normalisation is arbitrary. The DBI-Starobinsky 
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10" 



eain ^ 2.5 x 10-^ Sun ^ l- x W\ 7;^ ~ 5, 0i„/Mpi ~ 77.4193, Hin/A/pi 
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1.09 X 10^^ 

^■^ and j4_ = 0.92j4+. To make the comparison easier we have considered rig = 1 

0.96. Notice that working with non-negligible values of ei and/or £2 in the DBI 



and iVo = 10, where A'o is the number of e-folds at which the field goes through the transition. This implies A 

00/Mpi ~ 77.3989 and ^+/M|i ~ 2.368 x 10 ' 

for the Planck best fit while it is in fact n^ : 

Starobinsky model would lead to a significant tilt of the spectrum. 



amplitude of the oscillations would obviously be too large (compare for instance Figs. 6 and 9). In addition, 7 should 
not be too large. This has important implications for the exploration of the parameter space since it means that the 
relevant regime is in fact the "non-perturbative" one i.e. the one in which the analytical expressions derived above 
are no longer applicable. In this situation, only our mode by mode code can be used to explore the compatibility of 
the model with the data. 

Of course, another interesting aspect of the model is to which extent it produces non-Gaussianities. As is well- 
known, for slow-roll single field inflation with a standard kinetic term, the level of non-Gaussianity is very small, of 
the order of the slow-roll parameters, see Refs. [46-52]. In the case of the DBI-Starobinsky model, this is obviously 
no longer true. An interesting feature of the model is that, a priori, non-Gaussianities arise not only from one term, 
as is usually the case for non-slow-roll models, but from two (or even various) origins; the fact that the kinetic term is 
non standard and the discontinuity of the derivative of V{(j)) are the main sources of non-Gaussianity in this scenario. 
Concretely, the third order action reads [42, 53] 



M, 



PI 



drjdx 



2a 



eiSi 



ZHcl \ iclex 



/3 



'-lC5^P5^'id.,d,xdpd,x - 2a'-^C'5'^d^CdpX 
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a^e^ 



^ + 3Jcc- + ^^f§)cV 
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(69) 



[aEC7(jy'M2j)], E = eiH^M^^/cl and cx 



with X — 9 ^ 

M = 1 — 1/c^. In slow-roll canonical inflation 



-X/H'^ {dH/dX). 
the first vertex is absent because u 



The parameter u is defined by 
= i5i = 0. In DBI inflation, this 
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is no longer the case and it gives rise to non-vanishing non-Gaussianities with /J^^ ~ 35u/108 [42, 53]. Usually, as 
already mentioned above, the other contributions are negligible. On the other hand, if the kinetic term is minimal 
and the potential derivative has a discontinuity, then the third vertex proportional to £2 is the dominant one. Thus 
an important new feature of the scenario studied in this paper is that these two vertices are present. Moreover, 
one can expect the third vertex to be enhanced by the factor l/c^ since Cg ^ 1. The same mechanism should be 
valid for the other vertices as well. We are therefore in a rather complicated set-up in which non-Gaussianities are 
not easy to calculate since all terms contribute and since the mode function has a complicated behaviour. Despite 
that, it seems clear that the /i^l's parameters will be quite large, especially compared to the Planck constraints: 
/^°= = 2.7 ± 5.8, f^l = -42 ± 75 and /°["'° = -25 ± 39 [2]. For instance, using the DBI equation, one obtains 7 < 12. 
But, clearly, in our case, the constraint should be much tighter since other terms will contribute. In addition our 7 
is a time-dependent quantity so the calculation of the contribution of the first term will be modified. To estimate 
quantitatively the value of /^^ is a question that should be addressed by means of numerical calculations, or maybe 
using the formalism recently developed in [54] . It does not come as a surprise since we have shown before that already 
the two-point correlation function is an object difficult to calculate. We conclude that non-Gaussianites will be a very 
important probe to constrain the DBI-Starobinsky model. 

To end this paper, let us indicate the main directions for future works. Based on the previous considerations, it 
seems clear that the most promising direction is the calculation of non-Gaussianities. On the theoretical side, we have 
a new situation, not envisaged before, in which not only one vertex contributes but many and in a "coupled fashion" , 
i.e. the fact that Cg ^ 1 enhancing the contribution coming from the discontinuity of V. From the observational 
point of view, given the Planck result, it is clear that the corresponding constraints on the parameters of the model 
will be very tight. On the other hand, we have seen that the amplitude of the superimposed oscillations can be large 
(or, at least, seems larger than in the CS-model), even if the parameters are relatively close to standard slow-roll 
inflation. As a consequence, a priori, it remains possible that non-negligible superimposed oscillations improve the fit 
to the CMB data (especially by matching the Planck anomalies) while, at the same time, equilateral non-Gaussianities 
remain within the observational bounds. We hope to address this issue in more detail in the near future. 
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